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DETECTION OF PRESENCE OF CHEMICAL 
PRECURSORS 

ORIGIN OF THE INVENTION 

This application is a continuation-in-part of a prior-filed 
application, U.S. Ser. No. 11/178,079, filed 8 Jul. 2005. 

FIELD OF THE INVENTION 

This invention relates to use of functionalized carbon nano- 
structures to detect presence of one or more chemical precur- 
sors for a target molecule. 

BACKGROUND OF THE INVENTION 

Certain selected chemicals associated with terrorist activi- 
ties are too unstable to be prepared in a final form. These 
selected chemicals are often prepared as precursor compo- 
nents, to be combined at a time immediately preceding a time 
of application of the selected chemical. An example is a liquid 
explosive, which usually requires provision of an oxidizer, an 
energy source and a chemical or physical mechanism to com- 
bine the other components at a time immediately preceding 
detonation. Detection of presence of the oxidizer (e.g., H 2 0 2 ) 
or the energy source (e.g., nitromethane) is often possible but 
must be performed in a short time interval (e.g., 5-15 sec) and 
in an environment with a very small concentration (e.g., 
1-100 ppm), because the target chemical(s) is present in a 
sealed container. 

What is needed is a system that allows detection of pres- 
ence of a target oxidizer and/or a target energizer in small 
concentrations (as small as 1 ppm) in a relatively small time 
interval, preferably no more than about 5-15 sec. Preferably, 
the system should allow detection of at least one oxidizer and 
of at least one energizer, substantially simultaneously, should 
operate with a relatively small “footprint” in a real life envi- 
ronment, and should operate with only a small energy expen- 
diture. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
system and associated method for detecting one or more 
chemical precursors (components) of a multi -component 
compound that may become unstable when fully assembled 
or combined. First and second carbon nanostructures 
(“CNSs”) are loaded (by doping, impregnation, coating or 
other functionalization process) with different first and sec- 
ond chemical substances that react with first and second 
chemical precursors, respectively, which may be the same or 
may be different, if these precursors are present in a gas to 
which the CNSs are exposed. After exposure to the gas, a 
measured electrical parameter value EPV (e.g., impedance, 
conductivity, capacitance, inductance, etc.) changes with 
time in a predictable manner, if a selected chemical precursor 
is present, and will approach an asymptotic value promptly 
after exposure to the precursor. The measured EPVs are com- 
pared with one or more sequences of reference EPVs for one 
or more known target precursor molecules, and a most prob- 
able concentration value is estimated for each of one, two or 
more target molecules. An error value is computed, based on 
differences for the measured and reference EPVs using the 
most probable concentration values. Where the error value is 
less than an error value threshold, the system concludes that 
the target molecule is likely. Presence of one, two or more 
target molecules in the gas can be sensed from a single set of 
measurements. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates a system for practicing the 
invention. 

5 FIG. 2 graphically illustrates different concentration 
regimes. 

FIGS. 3 and 4 are flow chart of procedures for practicing 
the invention. 

FIGS. 5A, 5B and 5C graphically illustrate normalized 
10 responses to selected CNS loadings, before and after expo- 
sure to H 2 0 2 , H 2 0 and CH 3 OH. 

FIGS. 6A, 6B and 6C are histograms of response changes 
for the respective CNS loadings in FIGS. 5A, 5B and 5C, 

15 DESCRIPTION OF BEST MODE OF THE 

INVENTION 

FIG. 1 schematically illustrates a system 11 for practicing 
2Q the invention. The system 11 includes: a substrate 12; an 
appropriate catalyst underlayer 13; a sequence of sub arrays 
14-« of carbon nanostructures (“CNSs”) grown on the cata- 
lyzed substrate; a selected loading (doping, impregnation, 
coating, etc.) 15-72 of the CNS sub-array 14-72; a source 16 of 
25 a gas to be interrogated; a measurement mechanism 17 for 
measuring an electrical parameter value EPV (n) (impedance, 
conductance, capacitance, etc.) at each of the sequence of 
CNS sub-arrays, before and after exposure of that sub-array 
to the gas; and a computer 18 programmed to receive the 
30 sequence of measured electrical parameter values EPV (n), to 
compare the measured EPV(n) with a corresponding refer- 
ence values EPV(n;ref), and to estimate whether a target 
molecule is likely present in the gas in a most probable con- 
centration value. Each CNS sub-array 14-72 may consist of a 
35 single CNS or may include two or more CNSs, and two 
different CNS sub-arrays, 14-/2 and 14-72*, may have the same 
number, or a different number, of CNSs. The particular elec- 
trical parameter value EPV measured for each of the func- 
tionalized sensors may be electrical impedance, electrical 
40 conductance, capacitance, inductance or some other relevant, 
measurable electrical value. 

For relatively large concentrations, it is assumed that the 
measured EPV for a given target molecule will vary linearly 
with the concentration C, 

45 EPV 0 Q=a+b-C ; (1) 

where the parameters a and b are characteristic of the particu- 
lar target molecule for which the measurements are made. 
Where EPV(C lsl ) and EPV(C 12 ) are values for a pure sub- 
50 stance measured at the different concentration values C 1 1 and 
Ci 2 , respectively, the parameter values a and b can be esti- 
mated by 

a=(EPV(C 14 )C 1>2 -EPV(C 1 JC 14 }/{C lf2 -C 14 }, (2-1) 

55 

6=(EPV(C 14 )-EPV(C 1?2 }/{C 12 -C 14 }, (2-2) 

for each target molecule. For smaller concentrations («5-50 

ppm), it may be preferable to use a logarithmic approxima- 

„ tion, 

60 5 

EPV(C)=«'+^ , -log e C J (3) 

where a* and b* are selected parameters. The parameter values 
a* and b* can be determined in a manner similar to that of Eqs. 
65 (2-1) and (2-2), by replacing the variables C M and C 1#2 by 
log e {C M } and log e {C 12 }, respectively. FIG. 2 graphically 
illustrates the two concentration regimes for EPV(C), corre- 
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sponding to Eqs. (1) and (3), which join together at a concen- 
tration transition value C tr for which 

a'+b'-iog^C^a+b-C^. (4) 

Given an array of N sub-arrays of CNS sensors and a set of 5 
K target molecules (k=l, . . . , K;K^2), the EPV data for the 
sensors can be pre-processed in order to identify more clearly 
which CNS sensor sub-arrays are more sensitive to presence 
of a particular target molecule. In a first embodiment, the 
array is exposed to a selected target molecule, such as H 2 0 2 , io 
at a selected sequence {C } (q=l, . . . , Q with Q=N) known 
(not necessarily distinct) concentration values (e.g., 500 ppm, 
14,000 ppm, 65 ppm, 1200 ppm, 10 ppm, etc.), and a refer- 
ence measurement value EPV(n;q;k;ref) (n 1, 2, . . . , N) is 
recorded for each sensor sub-array n, each concentration C q 15 
of a reference gas containing a selected target molecule (k). 
Measurements for one or more concentration values for the 
selected target molecule gas can be repeated, if desired, to 
obtain an NxN square matrix of values. With q fixed, an NxN 
matrix {EPV(n;q;k;ref)/EPV(n;k;norm)}=E(n;q;k;ref) is 2 o 
formed, where, EPV(n;k;norm) is a normalization factor for 
the selected target molecule gas no. k and the sub-array no. n, 
which may be chosen as 


4 

event (degeneracy) that two or more of the N eigenvalues X(k) 
are equal, different non-zero linear combinations of the cor- 
responding eigenvectors V (k;X(k)) can be constructed that are 
orthogonal to each other, within a sub- space spanned by the 
reduced set of eigenvectors corresponding to the identical 
eigenvalues. 

Each matrix equation (9-k) has a sequence of N (eigen- 
value; eigenvector) pairs, {(X„(k);V„(k;X (k))} M , for a fixed 
reference gas k, and it is assumed here that the eigenvalues are 
arranged so that 

IX#)li=IX 2 (£)l^ . . . ^IXjv(A)I (10-k) 

and so that the highest magnitude eigenvalue in each set 
satisfies 

IX 1 (*1)I^IX 1 (*2)I^ . . . Ikfktyl, (11) 

where {kl, k2, . . . , kN} includes each of the integers 
{1, 2, . . . , N} precisely once. 

The eigenvector V(kl;X 1 (kl)) is identified as a first basis 
vector V'(l): 

V(k\;\(k\))=V(k\;\(k\)) (12-1) 


25 

(n y/p (5) 

EPV{n\ k\ norm ) = < ^ u q EPV{n\ q\ k\ ref) p > , 


where {u } is a selected set of non-negative weight numbers 
whose sum is a selected positive number (e.g., 1 or N) and p 
is a selected positive number. The normalization factor E(n; 
k;norm) may be a single EPV (e.g., u ql =l and u^O for q^ql), 
or may be a weighted sum of two or more EPVs. This nor- 
malization (optional) is intended to compensate for the con- 
centration dependence of the particular target molecule. 

For each sub-array n (fixed), the mean and standard devia- 
tion for each of K reference gases, numbered k=l, . . . , K 
(Ki^2) are computed, as 


N 

p{n\ k) = ^ E(n; q\ k; ref)/ N, 

< 7=1 

K 

Pin) = k)/K, 

k=l 

o-{n\ k) = I 'Yj re f) - *)} 2 in 
U= i 

N 

cr(n) = ^ cr(n; k) / K. 

< 7=1 

One now forms K auto scaled NxN matrices, defined by 

S(n;q;k)={E(n;q;k;ref)-[i(n) }/ o(n), 

and analyzes K eigenvalue equations 

S(n;q;k) V(k;k(k))=k(k) V(k;k(k)), 

where k (=1, . . . , K) is fixed and V(k;X(k)) is a normalized 
Nxl vector that will usually depend upon the reference gas 
(k) . If, as is likely, the N eigenvalues X(k) for a fixed reference 65 
gas k are distinct, the corresponding eigenvectors V(k;X(k)) 
are mutually orthogonal (non-degeneracy). In the unusual 


(8-k) 

60 

(9-k) 


40 

(6A) 

45 

(65) 

y 2 a a) 

I ■ 

(7B) 

55 


A second modified vector 

F(A:2;X 1 (A:2))=r(k2;X 1 (A:2))-{r(k2;X 1 (k2)),r(kl;X 1 

(kl))}V(kVMkl)) (12-2) 

is computed, where {V(k2;X 2 (k2)),V , (l)} is the scalar prod- 
uct (also referred to as the inner product) of the vectors 
V(k2;A, 2 (k2)) and V'(kl). More generally, a pth modified vec- 
tor 


V'{kp\\i{kp)) = 


(12 -p) 


p - i 


V(kp-,Xim)~Yj iV«p-,\i(kp)), V’(kr,\dkr)))V\kr, X^kr)) 

r= 1 


is computed for p=2, . . . , K. The set of vectors {V^kr;^ 
(kr))}^ (r=l, . . . , K) is mutually orthogonal, in the sense that 
the scalar products satisfy 


{V’{kr\ Ai (kr)), V'{ks ; Ai (ks))} = 6 ns . (13) 

> 0(r = s) 

= 0 (r*s). 


Each of the set of vectors {V^kr;^! (kr))} r (i^l , . . . , K) is 
maximally independent of each of the other vectors in the set, 
in the sense of mutual orthonormality (Eq. (13)). Each vector 
V(kr;>q(kr)) will have relatively large (primary) contribu- 
tions from some of the sensor sub -arrays and will have 
smaller (secondary) contributions from the remainder of the 
N sub-arrays. The vectors V^kr^^kr)) identify a maximally 
independent set of linear combinations of EPV responses 
from the N sub-arrays that can be used to distinguish presence 
of one reference gas (target molecule kr) from presence of 
another reference gas (target molecule ks). For example, if the 
set of reference gases are H 2 0 2 , H 2 0 and CH 3 OH, N=3 and 
three matrix eigenvalue equations are to be solved in Eqs. 
(9-k) (k=l, 2, 3). More generally, presence or absence of any 
of K target molecules (K^2) may be estimated. 
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The linear combinations LC(kp) of EPV measurements for 
the different sensor sub-arrays correspond to modified prin- 
cipal components for the particular reference gases chosen. 
Choice of another set of another set of reference gases will 
result in a different set of modified principal components, 
although change of one or more concentration values within 
a reference gas may have little or no effect on the modified 
principal components. 

FIG. 3 is a flow chart of a procedure for practicing the first 
embodiment. In step 31, N sub-arrays (Ni^2) of carbon nano- 
structures (“CNSs”), numbered n=l, . . . , N, are provided on 
a substrate, where each CNS sub-array is doped, coated, 
impregnated or otherwise functionalized (“loaded”) with a 
selected sensitizing substance drawn from a group of sub- 
stances, for example, Au nanoparticles, associated with other 
particles. In step 32, the CNS sub-arrays are exposed to a test 
gas to be interrogated (for presence or absence of K different 
target molecules, with KiU). In step 33, reference electrical 
parameter values E(n;q;k;ref)=(EPV(n;q;k;ref)/EPV(n;k; 
norm) are provided for each CNS sub -array no. n, for each of 
a selected sequence of K target molecules (KiU), and for 
each of a selected sequence of concentrations Cq of the target 
molecules (q=l, . . . , N). Here, EPV(n;k;norm) is a normal- 
ization factor that compensates for dependence of an EPV 
upon concentration of a target molecule. 

In step 34, means p(n) and standard deviations a(n) for the 
reference values E(n;q;k), summed over the concentration 
index q and over the target molecule index k, are provided. In 
step 35, autoscaled quantities 

S(n; q;k) ={E(n; q; k; ref)- \i(n ) }/ (a(« ) , (8-k) 

are provided. In step 36, (eigenvalue; eigenvector) pairs (X(k); 
V(k;X(k)) are provided that satisfy the matrix equations 

S{n;q;k) V(k;X(k))=X(k) V(k;\(k)). (9-k) 

The eigenvalues are assumed to be (re)arranged so that 

|X 1 (£)|^|X 2 (£)l^ ■ ■ ■ ^IVk)l(k=l, . . . , K) (10-k) 

and so that the highest magnitude eigenvalue in each set 
satisfies 

IX 1 (kll^lX 1 (k2)l^ . . . l?4(ftV)l, (11) 

where {kl, k2, . . . , kN} includes each of the integers 
{1, 2, . . . , N} precisely once. 

In step 37, K linear combinations V'(kp;X(kp))=LC{V(kr; 
X(kr))} =LC(kp) of the eigenvectors are formed so that a 
selected eigenvectors V’(kr;X 1 (kr)) are orthogonal to each 
other, in the sense of Eq. (13). This step may be implemented, 
for example, using an orthogonolization process set forth in 
Eqs. (12-1) through (12-p), or by any other suitable process. 
In step 38, the vectors {V , (kp;X 1 (kp))} /7 are interpreted as a 
maximally independent set of linear combinations of the 
EPVs for the N sub-arrays that distinguish between the K 
target molecules. Steps 39-44 (optional) enhance this 
embodiment by estimating whether a target molecule is likely 
to be present in a test gas. 

In step 39, normalized test measurements EPV(n;q;k;test)/ 
EPV(n;k;norm)=E(n;q;k;test) are provided for the (un- 
known) gas of interest. In step 40, linear combinations LC(kp; 
test) of these measurements corresponding to the linear 
combinations LC(kp;ref) of the N sub-arrays for the vectors 
{V^kpA^kp))}^ are provided. In step 41, error values 

e(kp)=I.e n \E{n;q;kp;test)-E{n;q;kp;ref) T, 
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are computed, where “LC(kp;ref)” indicates a weighted sum 
(w w ) over the linear combination of sensor sub-array indices 
(n) corresponding to the vector combination V^kp^^kp)) in 
Eq. (12-p) (p=l, . . . , K) and s is a selected positive number. 
5 In step 42, the error value e(kp) is compared with a selected 
positive error threshold value €(kp;thr), for at least one value 
of p=l, . . . , K. When €(kp) is less than €(kp;thr), the system 
interprets this condition as indicating that a target molecule, 
corresponding to the linear combination LC(kp) of sub-array 
to measurement values is likely to be present in the test gas, in 
step 43. When €(kp) is at least equal to €(kp;thr), the system 
interprets this condition as indicating that a target molecule, 
corresponding to the linear combination LC(kp) of sub-array 
measurement values is not likely to be present in the test gas, 
15 in step 44. 

In a second embodiment, illustrated for convenience here 
for first and second target molecules (K=2), each of the N 
CNS sensor sub-arrays is loaded (doped, impregnated, 
coated, or otherwise functionalized) with a different fimction- 
20 alizing substance FS, with each FS being chosen so that, for 
different first and second target molecules, at least one of the 
CNS sensor sub-arrays CNS no. n will produce substantially 
different EPV measurements for the first and second target 
molecules. Where conductance is the EPV of interest, the 
25 CNS loading substances FS may, for example, be Au nano- 
particles with associated side groups. 

Most probable concentrations, Cl (opt) and C2(opt), for the 
first and second target molecules present are computed using 
the approximations for concentration dependence set forth in 
30 Eqs. (1) and (3). For concentrations Cl and C2 above the 
transition values, Cl tr and C2 tr , a compound error value 


S(CU C2) = (15) 

35 

N 

Y K\ Epy ( n ) - R(l) + Ml) ■ Cl} - H( 2) + M2) ■ C2}| p , 


40 is computed, where {a„(l),b w (l) and (a w (2),b w (2)} are the 
parameter pairs (assumed known through experiment) that 
appear in Eq. (1) for the concentration dependence of the 
EPVs for the CNS sub-array no. n (n=l, . . . , N), w' w is a 
weighting number for the CNS sub-array number n, and p is 
45 a selected positive number. The weighting numbers w'„ may 
be equal or may reflect the relative number of CNSs in each of 
the sub-arrays. Compound error value in Eq. (15) is mini- 
mized with respect to choices of the (unknown) concentra- 
tions Cl and C2: 

50 

d{e(Cl;C2)}/dCl=0 at (Cl,C2)=(Cl(opt),C2(opt)) (16-1) 

d{e(Cl;C2)}/dC2=0, at (Cl,C2)=(Cl(opt),C2(opt)) (16-2) 

Using the linear concentration dependence from Eq. (1 ) for 
each of the target molecules 1 and 2, Eqs. (16-1) and (16-2) 
can be re-expressed as a first set of coupled relations in the 
optimum values, Cl=Cl(opt) and C2=C2(opt), 


0U N N N (17-1) 

Y W 'n b 2 n Cl + Y W n b n b ’n C1 = Y W n{ EPV i n ) ~ a n ~ <}M 

12=1 22= 1 22=1 

N N N (17-2) 

Y + Y W n b n C1 = Yj W J EPV W ~ a n~ a 'nK- 

65 22=1 22=1 22=1 


LC(kp;ref) 


(14) 
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Using the linear concentration dependence from Eq. (1 ) for 
target molecule 1 and the logarithmic concentration depen- 
dence from Eq. (3) for the target molecule 2, Eqs. (16-1) and 
(16-2) can be re-expressed as a third set of coupled relations 
in the optimum values, Cl=Cl(opt) and C2=C2(opt), 

N N N (17-3) 

X w' n blCl + X w’ n b n b'„\og'Cl = Y, <{EPV(n) -a,- a'„)b n , 

22=1 22=1 22= 1 

N N N (17-4) 

X < b nKCl + X < b >SeC2 = X <{EPV(n) -On-a'M, 

22=1 22= 1 n = 1 


where a common factor of C2 for all terms in Eq. (8-4) has 
been canceled. 

Using the logarithmic concentration dependence from Eq. 

(3) for the target molecules 1 and 2, Eqs. (16-1) and (16-2) can 
be re-expressed as a second set of coupled relations in the 20 
optimum values, Cl=Cl(opt) and C2 C2(opt), 

N N (17-5) 

^ H^^logpCl w^bnb^log^d = 25 

n=l n = 1 

N 

^ W n {EPV(n) -a n - a’ n }b n , 

n—l 

N N (17-6) 30 

w^bnb^lo^Cl + w^b'^log^d = 

n = 1 n = 1 

N 

^ W n {EPV{n) -a n - a' n }b ' n , 


The equation pair (17-1) and (17-2) is formally similar to the 
equation pair (17-5) and (17-6) so that the same well known 
algebraic methods can be used to formally determine the 40 
respective solution pairs (C1,C2) and (log e Cl,log e C2). 

Equations (17-1) through (17-6) extend, in an obvious 
manner, to Ki^3 different target molecules, where the EPV of 
each target molecule and each CNS sub-array (n) is repre- 
sented, in the appropriate concentration range, by Eq. (1) or 
Eq. (3). Where K(^2) target molecules are believed to be 
present, K coupled linear equations in the variable C or the 
variable log e C for these K target molecules are obtained, and 
the solutions provide most probable values Ck(opt) 5Q 
(k=l . . . , K) for the concentrations Ck in the gas. 

Preferably, Eqs. (17-1) and (17-2) are initially used to 
estimate the optimum values Cl (opt) and C2(opt). If the 
initial estimate of Cl (opt) is found to be at least equal to the 
first transition value Cl /r and the estimate of C2(opt) is found 55 
to lie below the second transition value C2, r in FIG. 2, Eqs. 
(17-3) and (17-4) would be used to redetermine the estimates 
Cl (opt) and C2(opt). If the initial estimates for Cl (opt) and 
for C2(opt) are found to lie below the respective transition 
values, Cl ?r andC2, r , Eqs. (17-5) and (17-6) would be used to 60 
redetermine the estimates Cl (opt) and C2(opt). If the initial 
estimates of Cl(opt) and C2(opt) are found to be at least equal 
to the respective first and second transition values Cl, r and 
C2, r , the initial estimates from Eqs. (17-1) and (17-2) would 
be used, unchanged, for the estimates Cl(opt) and C2(opt). 65 

The optimum value estimates, Cl (opt) and C2(opt), how- 
ever determined, are used to compute error values 


N 

(18-1) 

*1; Cliopt)) = X w„(1)|EPV(b) - EPV(Cl(opt); ref IF 1 , 
22=1 


N 

(18-2) 

e(2; C2(opt)) = X w„(2)| EPV{n) - EPV(Cl(optf ref2)\”\ 



22=1 


where w„(l) and w„(2) are non-negative weighting values, 
EPV(n) are the measured EPVs for the N-sensor sub -arrays, 
EPV(Cl) and EPV(C2) are reference EPVs, adjusted for con- 
centrations Cl and C2, for the first and second target mol- 
ecules, and pi and p2 are selected positive real numbers. 

The optimized error values, €(l;Cl(opt)) and €(2;C2(opt)), 
are then compared with selected error threshold values, e(l; 
thr) and €(2;thr), respectively, to determine if the target mol- 
ecule 1 and/or the target molecule 2 is likely to be present in 
the gas at the concentration Cl=Cl(opt) and/or at the concen- 
tration C2=C2(opt). The formalism set forth in the preceding 
is extendible to any number MiU of different target mol- 
ecules. At this point, the most probable concentration values, 
Cl (opt) and C2(opt), are determined, but a probability that 
the target molecules are present with these values is not yet 
estimated. 

FIG. 4 is a flow chart of a procedure for practicing the 
invention. In step 51, N sub-arrays (Ni^2) of carbon nano- 
structures (“CNSs”) are provided on a substrate, where each 
CNS sub-array is doped, impregnated, coated or otherwise 
functionalized (“loaded”) with a selected sensitizing sub- 
stance drawn from a group of substances, for example, Au 
nanoparticles, associated with other particles. In step 52, the 
CNS sub-arrays are exposed to a gas to be interrogated (for 
presence of K different target molecules, with KiU). In step 
53, a selected electrical parameter value EPV(n) for the CNS 
sub-array no. n is measured in a time interval of length At«5- 
30 sec or longer if desired . Optionally, the values EPV(n) for 
all sub-arrays are measured substantially simultaneously. 

In step 54, where K=2, a first reference value sequence 
EPV(n;refl) and a second reference value sequence EPV(n; 
ref2) are provided for sub-array no. n, corresponding to pres- 
ence of first and second different target molecules, present in 
a gas with known first and second concentration values, 
Cl(refl) and C2(ref2), respectively. 

In step 55, most likely concentrations, Cl(opt) and 
C2(opt), for the first and second target molecules in the test 
gas are estimated, based upon (1) differences between the N 
measured values EPV(n) of the collection of CNS sub-arrays 
and the corresponding reference values EPV(n;Cl;refl) and 
(2) differences between the N measured values EPV(n) of the 
collection of CNS sub -arrays and the corresponding refer- 
ence values EPV(n;C2;ref2), where Cl and C2 are unknown. 

At this point, the most-probable concentration values, 
Cl (opt) and C2(opt), for the first and second target molecules 
are known. In a first embodiment, it is sufficient to estimate 
the most-probable concentration values for the test gas. 

In an extension, an error value is computed, based on 
differences between the measured values EPV(n) for the N 
sub-arrays and computed concentration-dependent reference 
values, EPV(n;Cl«Cl (opt);ref 1) and EPV(n;C2«C2(opt); 
ref2), using Eq. (1) or Eq. (3) for the EPV for the first or 
second target molecule. Only if the first error value (or the 
second error value) is less than a selected first error threshold 
value (or less than a selected second error threshold value) is 
the most-probable concentration value Cl=Cl(opt) (or 
C2=C2(opt)) accepted as the concentration value for the first 
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target molecule (or the second target molecule) in the gas. In 
this extension, the procedures of steps 56-62 are optionally 
combined with the procedures of steps 51-55. 

In step 56, a first error value e(l) and a second error value 
€(2) are computed, depending upon (1) differences between 5 
the N measured values EPV(n) of the first CNS sub -array and 
the corresponding most -probable values EPV (n;Cl«Cl (opt) ; 
refl) and (2) differences between the N measured values 
EPV(n) of the second CNS sub-array and the corresponding 
most-probable values EPV(n;C2«C2(opt);ref2). 

In step 57, the first error value is compared with a first error 
threshold value €(l;thr). In step 58, when the first error value 
is less than the first error threshold value, the system interprets 
this condition as indicating that the first target molecule is 15 
present in the test gas with a selected first concentration 
Cl«Cl(opt). In step 59, when the first error value is at least 
equal to the first error threshold value, the system interprets 
this condition as indicating that the first target molecule is not 
present in the gas, or is present in the test gas with concen- 20 
tration substantially different from the most-probable con- 
centration, Cl«Cl(opt). 

In step 60, the second error value is compared with a 
second error threshold value e(2;thr). In step 61, when the 2 5 
second error value is less than the second error threshold 
value, the system interprets this condition as indicating that 
the second target molecule is present in the test gas with a 
selected second concentration C2«C2(opt). In step 62, when 
the second error value is at least equal to the second error 30 
threshold value, the system interprets this condition as indi- 
cating that the second target molecule is not present in the test 
gas, or is present in the gas with concentration substantially 
different from the most probable concentration, C2«C2(opt). 

Optionally, the system can test for presence of one of K^2 
different target molecules substantially simultaneously. One 
advantage of the invention is its flexibility: presence of any 
reasonable number of target molecules can be tested for with 
a single set of EPV measurements. 

FIGS. 5A, 5B and 5C graphically illustrate measurements 
of normalized (dimensionless) responses, 

NR(t)={R(t)-R 0 }/R 0 , (19) 

of electrical resistance of a loaded CNT, measured before and 45 
after introduction of a reference gas H 2 0 2 (5 A), H 2 0 (5B) and 
CH 3 OH (5C), into a chamber containing the loaded CNT, 
minus the measured resistance R 0 of an unloaded CNT, 
divided by R 0 . The curve marked “gas,” which rises substan- 
tially vertically beginning at a time t«0.43 min, indicates 50 
introduction of the reference gas into the chamber. As an 
example, the changes in NR(t) for C 12 , measured at times 
before and after introduction of the reference gas, are mea- 
surable: ANR(t)~- 0.05 (H 2 0 2 ), «-0.02 (H 2 0) and «+0.01 
(CH3OH) for the reference gases examined. 

FIGS. 6A, 6B and 6C graphically illustrate histograms 
representing re-normalized responses for the respective ref- 
erence gases H 2 0 2 , H 2 0 and CH 3 OH, computed as 

NNR{NR(t:before)-NR(t;after)}/NR(t;before), (20) 60 

for 8 sensors among 24 sensors in an embodiment (FIG. 1), 
where “before” and “after” have the meanings discussed in 
connection with FIGS. 6A, 6B and 6C. Note that most of the 
histogram amplitudes are positive, with the exception of the 65 
renormalized responses for sensors no. 2, 3 and 4 for H 2 0 and 
for CH3OH. 


What is claimed is: 

1 . A method for estimating presence of a target molecule in 
a gas, the method comprising: 

providing N carbon nanostructure (“CNS”) sub -arrays, 
numbered n=l , . . . , N (N^2) on a substrate, where each 
CNS sub-array is loaded with a selected sensitizing sub- 
stance comprising Au nanoparticles; 
estimating a first functional relationship EPV(C1;1) 
between a concentration Cl of a first target molecule and 
a first measured electrical parameter value EPV(n=l) in 
a gas, where the first relationship has a first approximate 
form EPV(n= l)=ar+bl' log(Cl) in a first concentration 
range of Cl and has a second approximate form EPV 
(n=l)=al+bl-Cl in a second concentration range of Cl, 
where al', bl', al and bl are selected real numbers; 
estimating a second functional relationship EPV(C2;2) 
between a concentration C2 of a second target molecule 
and a second measured electrical parameter value EPV 
(n=2) in a gas, where the second relationship has a third 
approximate form EPV(n=2)=a2'+b2 , *log(C2) in a first 
concentration range of C2 and has a fourth approximate 
form EPV(n=2)=a2+b2*C2 in a second concentration 
range of C2, where a2', b2', a2 and b2 are selected real 
numbers; 

exposing the CNS sub -array no. n (n=l, . . . , N) to the gas, 
and measuring the electrical parameter value EPV(n) in 
a time interval of length no more than about 15 sec; 
estimating a probable value of concentration, Cl(opt) and 
C2(opt), of the first target molecule and the second target 
molecule, respectively, in the gas, based upon (1) differ- 
ences between the value EPV(n=l) and the first func- 
tional relationship for the first target molecule and upon 
(2) differences between the value EPV(n=2) and the 
second functional relationship for the second target mol- 
ecule 

computing a first error value depending upon differences 
40 between the measured values EPV(n) of the CNS sub- 
array no, n=l and reference values E(Cl(opt);ref;n=l) 
that corresponds to presence of the first target molecule 
in a gas with a concentration of Cl (opt); 
computing a second error value depending upon differ- 
ences between the measured values EPV(n) of the CNS 
sub-array no, n=2 and reference values E(C2(opt);ref; 
n=2) that corresponds to presence of a second target 
molecule in a gas with a concentration of C2(opt); 
when the first error value is less than a selected first error 
threshold value, interpreting this condition as indicating 
that the first target molecule is present in the gas with a 
concentration of about Cl (opt); 
when the second error value is less than a selected second 
error threshold value, interpreting this condition as indi- 
cating that the second target molecule is present in the 
gas with a concentration of about C2(opt); 
when the first error value is at least equal to the first error 
threshold value, interpreting this condition as indicating 
that the first target molecule is not present in the gas or is 
present in the gas with a concentration substantially less 
than Cl (opt); and 

when the second error value is at least equal to the second 
error threshold value, interpreting this condition as indi- 
cating that the second target molecule is not present in 
the gas or is present in the gas with a concentration 
substantially less than C2(opt). 
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2. The method of claim 1, further comprising: 
computing said first error value as a sum 

N 

ei=Yj ™n(mPV(n) - E(Cl(opt); ref-, n = l)| pI , 

> 1=1 

where EPV(n) is said measured parameter value of said first 
CNS sub-array, EPV(Cl(opt);refl) is a corresponding refer- 
ence parameter value for said first target molecule, w„(l) is a 
non-negative weighting value for said first target molecule, 
and pi is a selected positive number. 

3. The method of claim 1, further comprising: 
computing said second error value as a sum 

N 

£2 = y w„(2)| EPV(n) - E(C2(opt)-, ref-, n = 2)\ p2 , 

n= 1 


12 

where EPV(n) is said measured parameter value of said sec- 
ond CNS sub-array, EPV(C2(opt);ref 2) is a corresponding 
reference parameter value for said second target molecule, 
w„(2) is a non-negative weighting value for said second target 
5 molecule, and p2 is a selected positive number. 

4 . The method of claim 1 , further comprising choosing said 
electrical parameter value EPV(n) from the group of param- 
eters consisting of electrical impedance, electrical conduc- 
tance, capacitance and inductance. 

5 . The method of claim 1 , further comprising estimating 
said first functional relationship where said probable value of 
concentration Cl (opt) is no greater than about 50 ppm. 

6 . The method of claim 5 , further comprising estimating 
1 5 said first functional relationship where said probable value of 

concentration Cl (opt) is no greater than about 5 ppm. 



